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An elementary approach for computing the values at nega- 
tive integers of the Riemann zeta function is presented. The 
approach is based on a new method for ordering the integers 
and a new method for summation of divergent series. We 
show that the values of the Riemann zeta function can be 
computed, without using the theory of analytic continuation 
and functions of complex variable. 
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1 Introduction 

The Riemann zeta function C{s) is one of the most important objects in the 
study of number theory (if not in the whole modern mathematics). It is a 
classical and well-known result that C{s)-, originally defined on the half plane 
3fJ(s) > 1, can be analytically continued to a meromorphic function on the 
entire complex plane with the only pole at s = 1, which is a simple pole with 
residue 1. 

One of the main reasons of interest to C^{s) is that the special values of 
(^(s) at integers have been proved or conjectured to have significant arithmetic 
meanings. For instance, Zagier's [IHl HI] conjecture concerning the relation 
between (^(n) and the n- logarithms forn > 2 and Lichtenbaum's [26] conjecture 
connecting (^(n) with motivic cohomology. 



* extended version of the contribution to the proceedings of the XXVII International 
Colloquium on Group Theoretical Methods in Physics, Yerevan, Armenia, August 13-19, 
2008 

'^Institution of the Russian Academy of Sciences, V.A. Trapeznikov Institute for Control 
Sciences of RAS, Moscow, Russia; E-mail: abagd asariOhotmail . com) 



1 



Although the Riemann zeta function is considered as being primarily rele- 
vant to the "purest" of mathematical disciplines, number theory, it also occurs, 
for instance, in applied statistics (Zipf's law, Zipf-Mandelbrot law), physics, 
cosmology [9l[T0l[TT]. For example, its special values C(3/2) is employed in cal- 
culating the critical temperature for a Bose- Einstein condensate, ({4) is used 
in Stefan-Boltzmann law and Wien approximation. The Riemann zeta func- 
tion appears in models of quantum chaos [H [5l |22l |3T] and shows up explicitly 
in the calculation of the Casimir effect [271 128] . 

The method of zeta function regularization [12], which is based on the 
analytic continuation of the zeta function in the complex plane, is used as one 
possible means of regularization of divergent series in quantum field theory 

So, the Riemann zeta function is an interesting object for study not only 
for mathematicians but for physicists as well. 

2 Zeta Function and its Values at Negative 
Integers 

The zeta function was first introduced by Euler and is defined by 



The series is convergent when s is a complex number with 3ft(s) > 1. 

In 1859 Riemann defined ({s) for all complex numbers s by analytic con- 
tinuation. Several techniques permit to extend the domain of definition of the 
zeta function (the continuation is independent of the technique used because 
of uniqueness of analytic continuation). One can, for example, consider the 
zeta alternating series, so-called the Dirichlet eta function 



defining an analytic function for 3?(s) > 0. When the complex number s 
satisfies 3f?(s) > 1, we have 





oo 



u=l 



oo ^ oo 2 2 

u=l u=l ^ ' 
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or 

C(^) = Y^^, ^is) > 1. 

Since r]{s) is defined for 9f?(s) > 0, tfiis identity permits to define the 
zeta function for all complex numbers s with positive real part, except for 
s = 1 + 2n7Ti/ log 2 where n is an integer. For n = 0, we have a pole at s = 1, 
and for n ^ one can use the derivative of the ?7-function since it is known 
that T] has zeros at these isolated points on the line 3fJ(s) = 1 [251 ESI ISU Bi] . 

So, the zeta alternating series is linked with the original series (the zeta 
function) by the simple relation 

(1) C{s) = {1 - 2'-^)C{s) 

where we have put ({s) = r]{s). 

Around 1740, Euler [131 [El [IS] discovered a method of calculating the 
values of the divergent series 



1+1+1+14 
1+2+3+4+ 



1 

12 

1 + 4 + 9 + 16 + ... = 

1 



+ 27 + 64 + ... 



120 
etc. 

In modern terms, these are the values at non-positive integer arguments 
of the Riemann zeta function C{s), which, as it was said above, defined by the 
series, absolutely convergent in 3fJ(s) > 1 

111 
C(^) = l + ^ + ^+4l + - 

To find these values, Euler considered the alternating series 
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(2) 1™- 2™ + 3" -4™ + 5™ -6" + ... 



He observed that the value ([2]) is obtained as a limit of the power series 
l'"-2'"x + 3'"x2-4'"x3 + ... 
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as X — )■ 1, since, although the series itself converges only for |x| < 1, it has 
an expression as a rational function (analytic continuation, as we now put it), 
finite at x = 1, which is obtained by a successive application of multiplication 
by X and differentiation (equivalently, applying the Euler operator x ■ d/dx 
successively after once multiplied by x) to the geometric series expansion 

(3) Y^—^ = l-x + x'^-x^ + x*-x^ + ... (|a;| < 1) 

For instance, if we substitute a; = 1 in ([3]), we find formally 

l = l_l + l_l + l_l + ... = ^(0) 

and hence, in view of ([T]), we have C(0) = — |. More examples are 

^ = 1 - 2x + -4a;^ + 5a;^ - ... 

1 - 2^^ + 3^x2 - A^x^ + 5V - ... 



(l + x)2 

1 — X 



'1 + x 



l3 



1 — 4x + x^ 
il + xY 



2^x + 3^x2 - 4^x^ + 5^x^ - 



which give us 



and in turn 



C(-l) = ^, C(-2) = 0, C(-3) 



C(-l) = -i^, C(-2) = 0, C(-3) 



12 ' ' 120 

Euler found these values of divergent series, having no notion of the ana- 
lytic continuation and functions of complex variable. 

In modern concepts this method provides no rigorous way for obtaining 
the values of at non-positive integers because it is commonly considered 
that the values of (^{—m) should be established as values of the analytically 
continued function C,{s) at s = — m. 

There are different approaches to evaluation of the values of C,{—m). These 
values can be computed by using the functional equation satisfied by this 
function [33], by using the Euler- Maclaurin formula [ISIE31, in [3D] the values 
of C(— m) are derived from a particular series involving sC,{s + 1), another 
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method is presented in [29]. Other approaches are based on the so-called 
g-extensions (g-analogs) of Riemann zeta function [211 EHl El]. 

However, in this paper we find in a mathematically rigorous way a new ap- 
proach to the values of ({s) at non-positive integer points, without addressing 
the theory of analytic continuation and functions of complex variable. 



3 Definitions and Preliminaries 

We shall deal with two series 

= / (Dirichlet eta function) 

and 



u=l 



Cl-^) = / — (Riemann zeta function) 

u=l 

and aim to evaluate C^{s) and C,{s) at non-positive integers {C,{—m),C{—'m)) in 
elementary and real way, without using complex- analytical techniques. 

As a starting point we take a new method for ordering the integers [371 
[38] , which provides very well not only a real and elementary approach to 
computing the values of C(s) at negative integers but also a potentially new 
field of research [51 ESI EB] • 

In this paper we restrict ourselves to considering the values of C,{s) and 
C(s) at non-positive integers and leave the issue of computing the values of 
C(s) and C(s) at positive integers for our next work [7j- 

To make paper self-contained, we introduce some basic definitions and 
propositions [55] (proofs can be found therein), necessary for the aim of this 
paper. 

To be more specific, let us consider the set of all integer numbers z|3 

Definition 1 We say that a precedes b, aJ) G Z , and write a ~< b, if the 
inequality < holds; a -< b <^ ^ < 

This method of ordering, obviously, gives that any positive integer number 
(including zero) precedes any negative integer number, and the set Z has 
zero as the first element and —1 as the last element, i. e. we have Z = 



^we denote by a, b, c integer numbers and by n, m, k natural numbers 
^assuming by convention 0^^ = cxd 
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[0,1,2,... — 2,-1] §. In addition, the following two necessary conditions of 
axioms of order hold: 

1. either a -< b ot b -< a 

2. if a -< 6 and 6 -< c then a -< c 

Definition 2 A function f{x), x G Z , is called regular if there exists an 
elementary function F{x) such that F{z + 1) — F{z) = f{z), Wz G Z . The 
function F{x) is said to be a generating function for f{x). 

Remark 3 If F{x) is a generating function for f {x) , then the function F {x) + 
C, where C is a constant, is also a generating function for f{x). So, any 
function F{x) which is generating for f{x) can be represented in the form 
F{x) + C{x), where C{x) is a periodic function with the period 1. 

Suppose f{x) is a function of real variable defined on Z and Za,b is a part 
of Z such that Za^b = [o-, b] if a ^ b and 5 = Z \ {b,a) if a >- b, where 
Z\{b,a) = [a,-l] U [0,6]. 

Definition 4 For any a,b & Z 

(4) E/H= E /("fl- 

U=a 'fJ'&^a.b 

This definition satisfies the condition of generality and has a real sense 
for any integer values of a and b (a<6). The definition H] extends the classical 
definition of sum Ylt /(^) '^^^^ b < a. 

We introduce for regular functions the following quite natural conditions: 



n 00 [~| 

1. lfSn=E /(«) Vn, then lim Sn = E fim 



n/2 00 

2. If Sn= Yl /(«) then lim Sn = Yl /(^) 

u=l u=l 



3. If = '5', then Y (^f{u) = aS, a E R. 



•^geometrically, the set Z can be represented as cyclically closed (closed number line) 
"^observing the established order of elements Za^b 

— > 00 means that n unboundedly increases, without changing the sign 
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oo oo oo 

4. If E = S, and E f2{u) = S^, then E + /2(w)) = + ^2. 

u=a u=a u=a 



h 

5. For any a and b,a<b: F{h + 1) - F(a) = ^ f{u). 

u=a 

6. If G = [ai, hi] U [a2, ^2], [oi, hi] fl [02, ^2] = 0, then 

u£G ii=ai u=a2 

The conditions ([I])-© define a method of summation of infinite series, 
which is regular due to ([T]). 

Proposition 5 // f{x) is a regular function and a E Z is a fixed number, 
then 

a-l 

(5) E/H = E/(") 

u=a ueZ 

Proposition 6 For any numbers m and n such that m -< n 

n —m 

(6) E = E 

u=m u=~n 

Proposition 7 Let f{x) be a regular function and let a, h, c be any integer 
numbers such that h G Z^^c- Then 

c b c 

(7) E/H = E/H+ E' 

u=a u=a u=b-\-l 

Proposition 8 Suppose f{x) is a regular function. Then 

a-l 

(8) ^/(m) = VaGZ 

u=a 

or, which is the same in view of ^ 

(9) E/(^) = 

c 

^the prime on the summation sign means that ^ fi"^) = for = c 

u=b+l 
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Proposition 9 For any regular function f{x) 



a-l 



(10) ^/(n) = - J2 /H Va,6eZ 

u=a u=b+l 

From ( JTOl) . letting a = and 6 = —n, we have 

-n. -1 

«=0 «=— n+1 

and using ([6]), we get 

— n 71—1 

u=0 u=l 

and 

— n 71—1 

(11) E/(^) = -E/(-^) 

Using (|TT|) . we obtain 
Theorem 10 For any even regular function f{x) 

(12) £/(-)= 

u=l 

independently on whether the series is convergent or not in a usual sense. 

Example 11 Consider some examples of both convergent and divergent 
ries. 



1. Convergent series 

oo 

y ^ =1 

^4^2-1 2 



u=l 

(2^^-1/2)2 
^ (4"-1)(m-1/2)-1 _ 1 

u=l 



E°° (m^ + 1/4) tan(l/2) cosM — -usiiiM tan(l/2) 
(4m2 - 1)2 " 8 

u=l ^ ' 

with the generating functions, respectively 

-1 f-1)""^ 
F{n) = — -, F{n' ^ ' 



2(2n- 1)' ^ ' (2n- 1) 



2' 



p,.^^{n-m ^ sin(n - 1/2) 

2"^-"+i ' 8(2n-l)2cos(l/2) 



2. Divergent series 



1 

1-1 + 1-1 + ... = - 
2 



1 + 1 + 1 + 1 + ... = -^ 



l2fc + 22'= + 32^.. = VA; 

wi/i i/ie generating functions, respectively 

F(n) = ^^, F(n) = n-1, F(n) = 52fe(n - 

Theorem 12 For any polynomial f{x), x E R 
(13) lim (-1)"/H = 





(14) lim f{n) = [ f{x)dx 

-1 
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4 Main Results 

Using (flSll we immediately obtain 



Theorem 13 Let a{x) and be elementary functions defined on Z and 

satisfying the condition a{x) — f3{x) = f{x), where f{x) is a polynomial. Sup- 
pose that fi{x) is a function such that fj,{x) = a{x) if 2 \ x and fj,{x) = /3{x) 
if 2\ X. Then 

lim jjL{n) = - lim (a{n) + f3{n)). 

From Theorems [12] and [131 we get the following 

Theorem 14 Let = ai + (u — l)d, d > 0, is an arithmetic progression. 
Then 

5d — Qai 



1) XI 



12 

u=l 



2) 5^(-ir-v = ^ 

Example 15 

oo ^ 

^1 = 1 + 1 + 1 + ... = --, {d = 0) 

u=l 

oo ^ 

^^ = 1 + 2 + 3 + ... = --, {d = l) 

u=l 

oo ^ 

J](2m-1) = 1 + 3 + 5 + ... = -, (d = 2) 



etc. 



oo ^ 

= + + = (rf = o) 

u=l 

oo ^ 

5^(-ir~^«=l-2 + 3-4 + ... = -, (d=l) 

u=l 

oo 

^(_1)«-1(2m _ 1) = 1 _ 3 + 5 _ 7 + ... = 0, (fi = 2) 



u=l 

etc. 
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Let us now consider the Bernoulli polynomials. Bernoulli polynomials can 
be defined, in a simple way, by the symbolic equality Bn{t) = {B + 1)", where 
the right hand side members should be expanded by the binomial theorem, 
and then each power 5" should be replaced by and Bn are the Bernoulli 
numbers (analogously, the symbolic equality Bn = {B -\- 1)" can be written for 
the Bernoulli numbers). 

Bernoulli numbers play an important role in many topics of mathematics 
like analysis, number theory, differential topology, and in many other areas. 
These numbers were first introduced by Jacob Bernoulli (1654-1705) and ap- 
peared in Ars Conjectandi, his famous treatise published posthumously in 
1713, when he studied the sums of powers of consecutive integers l'^ + 2^^ + 
3^ + ... + n^. They often occur when expanding some simple functions in a 
power series. For instance, in the series 

1 o2 D o4 D o2A; r) 
cot X = -^X + -^X^ - ... - (-i f -Px^'^-l + ... < \x\ < TT 

X 2! 4! 2kl 

which appeared in astronomical works of J. Bernoulli. 
Let 

u=0 ^ ^ u=\ 

be the Bernoulli polynomial. 

We derive two well-known equalities, without invoking the complex-analytical 
notions. 

Since Bk{n) — Bi^{n — 1) = n'^, the function f[x) = x^ is regular. Then, 
in view of ([TTl). we have 



(15) Bu{-n) = {-lf-'Bk{n 

On the one hand, according to f|T5l) 



On the other hand, in view of ( fl4l) 



fc-l 



/I /i- -I- 1 \ 

Bu-i{x)dx = ——-y^{-iT^H ^ ]B^ = y^u^-' 
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Adding these two, we get 

oo 

(16) 5^«'"' = -^ V/t 

u=l 

Now taking the formula 
(2 



U=l ^ ^ ?i=l 



k 



which holds for any natural number n, and then passing to the limit (letting 
n — )■ oo) and taking into account ( !T3|) . we obtain 



(17) 5^(-i)"-1m'=~i = -^^Bk yk 



u=l 



So, we have obtained the equalities ffTB]) [35j and ffTTl) [20], and all the 
above propositions solely within the framework of real analysis. 

Combining f|T6|) . f|T7|) . Theorems [TO] and [T^ and putting k — 1 = m, we 
elementarily and immediately arrive at 

Proposition 16 For each non-negative integer m, we have 
(i) for Riemann zeta function 



I'" + 2"^ + 3" + 4™ + ... = C(-m) 



Bm+l 



m + 1 

(a) for Dirichlet eta function 

^ rym+l 1 

1*" - 2"^ + a*" - 4™ + ... = C{-m) = Bm 

m + 1 

Using f|T5|) . we can also show that all odd Bernoulli numbers Bk (except 
for Bi = 1/2) are equal to zero. Indeed, 

p-l)/2] . , . 

B,(n) -i-ir-' B,{-n) = ^ J2 (2I + 1 j = 



u=0 



whence we immediately obtain i?i = 1/2 and i?2«-i = 0, m = 1, 2, .... 
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5 Conclusion 



In conclusion, we would like to emphasize that, unlike of the many other 
methods that make use of either contour integration or analytic continuation 
of C{s), our approach, elaborated within the new theoretical direction [38], has 
allowed us to obtain the values of the Riemann zeta function (and also the 
Dirichlet eta function) at non-positive integers in a purely real way, without 
any notions of complex analysis and analytic continuation, and it follows from 
one single concept. It would be interesting to apply the similar reasonings and 
techniques, within the setting of [5B], to the some zeta related functions and 
other special functions. We plan to study this in our subsequent works. 
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